PERMUTATIONS WITHOUT LONG DECREASING 
SUBSEQUENCES AND RANDOM MATRICES 

PIOTR SNIADY 

Abstract. We study the shape of the Young diagram A associated via 
the Robinson-Schensted-Knuth algorithm to a random permutation in 
Sn such that the length of the longest decreasing subsequence is not big- 
ger than a fixed number d; in other words we study the restriction of the 
Plancherel measure to Young diagrams with at most d rows. We prove 
that in the limit n — > oo the rows of A behave like the eigenvalues of a 
certain random matrix (namely the traceless Gaussian Unitary Ensemble 
random matrix) with d rows and columns. In particular, the length of the 
longest increasing subsequence of such a random permutation behaves 
asymptotically like the largest eigenvalue of the corresponding random 
matrix. 



1. Introduction 

1.1. Formulation of the problem. Let an integer d > 1 be fixed. For 
any integer n > 1 we consider the set of the permutations n E Sn such 
that the length of the longest decreasing subsequence of n is not bigger 
than d; in other words it is the set of the permutations avoiding the pattern 
( d + 1 , d, . . . , 3, 2, 1 ) . Let be a random element of this set (probabilities 
of all elements are equal). In this article we are interested in the following 
problem: 

Problem 1. Let G Sn be a random permutation with the longest de- 
creasing subsequence of length at most d. What can we say about the as- 
ymptotic behavior of the length of the longest increasing subsequence of Un 
in the limit r\. ^ oo? 

Let \n = (An,i, . . . , An,d) be the (random) Young diagram associated 
via the Robinson-Schensted-Knuth algorithm to tt^ (notice that since the 
number of the rows of A^^ is equal to the length of the longest decreasing 
subsequence of Un, A^ has at most d rows). In other words, An is a random 
Young diagram with at most d rows, where the probability of the Young 
diagram A is proportional to (dim Pa)^, where dim Pa denotes the dimension 
of the corresponding irreducible representation of Sn; therefore, if we drop 
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the restriction on the number of the rows of the Young diagrams (which 
can be alternatively stated as d > n), then the distribution of 'Kn is the 
celebrated Plancherel measure. 

Since An,i is equal to the length of the longest increasing subsequence in 
Tin, Problem [T] is a special case of the following more general one: 

Problem 2. What can we say about the asymptotic behavior of the random 
variables (A^^ i , . . . , A^^ d) the limit n ^ oo? 

1.2. Case d = 2. The first non-trivial case d = 2 was considered by 
Deutsch, Hildebrand and Wilf [IDHW03 II. In this case the random variables 
An,i , Ati,2 are subject to a constraint A^j + An,2 — n- therefore it is enough 
to study the distribution of An,i. Deutsch, Hildebrand and Wilf proved that 

the distribution of (Anj — 5) converges to the distribution of the length 

of a random Gaussian vector in M^; in other words ^ (An,i — j)^ converges 
to the xl distribution with 3 degrees of freedom (a careful reader may no- 
tice that the authors of [|DHW03il use a non-standard definition of the 

distributions and therefore they claim that l^nj — itself converges 

to xi)- Their proof was based on an explicit calculation of the number of 
the permutations which correspond to a prescribed Young diagram with at 
most two rows. 



1 .3. Case d = oo. Another extreme of this problem is to consider d = cxd; 
in other words, not to impose any restrictions on the random permutations 
Tin. In this case the random Young diagram A^ is distributed according 
to the Plancherel measure. The authors of IIBDJ991 lOkoOOl IJohOlU proved 
that the joint distribution of longest rows of A (after appropriate rescaling) 
converges to the same distribution (called Tracy-Widom distribution) as the 
joint distribution of the biggest eigenvalues of a large random matrix from 
the Gaussian Unitary Ensemble. 

1.4. The main result: intermediate values of d. We equip the vector 
space of d X d Hermitian matrices with a Gaussian probability measure 
with a density 

with respect to the Lebesgue measure, where Zd is the normalizing con- 
stant. We say that a random matrix (Aij)i<tj<d distributed accordingly to 
this measure is a Gaussian Unitary Ensemble (GUE) random matrix. 

We call B = A — -^TrAa traceless Gaussian Unitary Ensemble ( GUEq) 
random matrix; it corresponds to the Gaussian probability measure on the 
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set of d X d Hermitian matrices with trace zero and the density 

with respect to the Lebesgue measure, where is the normalizing con- 
stant. 

The joint distribution of eigenvalues for GUE is well-known [|Meh91|| . 
which allows us to find the corresponding distribution for GUEq; namely, 
if Xi > • ■ ■ > are the eigenvalues of a GUEq random matrix, then their 
joint distribution is supported on the hyperplane xi + ■ ■ • + Xd = with the 
density 

(1) T^e-^^n^'^^-^^) 

with respect to the Lebesgue measure, where Cd is the normalization con- 
stant. 

Theorem 3 (Main theorem). Let the integer d > 1 be fixed; for each n > 1 
let An = (ArL,i , . . . , An,d) be, as in Section li.il a random Young diagram 
with TL boxes and with at most d rows. 

Then the joint distribution of the random variables ( \ — [^nx ~ 



i<i<d 

converges, asn go, to the joint distribution of the eigenvalues of a GUEq 
random matrix. 

We postpone its proof to Section [2l 

Corollary 4. Let d > 1 be fixed, and for each n > ^ let n^^ E Sn be a 
random permutation with the longest decreasing subsequence of length at 
most d. We denote by An,i the length of its longest increasing subsequence. 

Then the distribution of (An,i — converges to the distribution of the 
largest eigenvalue of the GUEq random matrix. 

It should be pointed out that the distibution of eigenvalues of a GUEq 
random matrix appears also in a related asymptotic problem HJohOlU of the 
distribution of the rows of a Young diagram associated (via RSK algorithm) 
to a random word consisting of n letters in an alphabet of d symbols in the 
limit of n — > cxD. 

1.5. Case d = 2 revisited. The set of 2 x 2 Hermitian matrices with 
trace zero can be viewed as a 3-dimensional Euclidean space with a scalar 
product (A,B) = TrAB. A GUEq random matrix can be viewed under 
this correspondence as a Gaussian random vector in the coordinates of 
which are independent with mean zero and variance 1 . Each 2x2 Her- 
mitian traceless matrix A has two eigenvalues xi = A, X2 = —A, where 
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^ ~ 72 1 1 1 1 ^ V "^^T^- Therefore, for a GUEo random matrix the corre- 
sponding random variable 2(xi )^ is distributed like the square of the length 
of a standard Gaussian random vector in M^, which coincides with the xi 
distribution; thus for d = 2 Corollary |4] allows us to recover the result of 
Deutsch, Hildebrand and Wilf [IDHW03L 



1.6. Idea of the proof . In Section|2]we will prove Theorem[3l the main re- 
sult of this article. Our proof will be based on an explicit calculation of the 
number of standard Young tableaux with a prescribed shape. The standard 
method to do this would be to use the hook-length formula, which would 
be not convenient for our purposes. Instead, we will use the determinantal 
formula of Frobenius and MacMahon. In order to make the connection to 
random matrices more explicit we shall recall its proof due to Zeilberger 
llZei83ll which is based on the observation that a Young tableaux with at 
most d rows can be viewed as a certain trajectory of d non-colliding parti- 
cles on a line. Thus we will find explicitly the asymptotic joint distribution 
of the rows of a Young diagram; this distribution turns out to coincide with 
the distribution ([T]) of the eigenvalues of a GUEq random matrix. 

The reader may wonder if the connection between Young diagrams and 
random matrices given by Theorem [3] might be purely accidental. In the 
following paragraph we will argue why it is not the case and how deep 
connections between Young diagrams and random matrices may be seen in 
our proof of Theorem |3l 

In the above discussion we treated the distribution ([T]) of the eigenvalues 
of a GUEo random matrix as granted; now let us think for a moment about 
its derivation. GUEq is a Gaussian matrix; for this reason (up to a simple 
scaling factor) it can be viewed as a value at some fixed time of a matrix- 
valued Brownian bridge. It is known [ Dys62[ lGra99l that the eigenvalues 
of a matrix- valued Brownion motion behave like Brownian motions condi- 
tioned not to collide. Since a matrix-valued Brownian bridge is a matrix- 
valued Brownian motion conditioned to be zero at time 1 , it follows that 
its eigenvalues form Brownian motions conditioned not to collide and to 
be zero at time 1 ; in other words these eigenvalues form Brownian bridges 
conditioned not to collide. In this way the determinantal formula of Karlin 
and McGregor [KM59| can be applied. In the conditioning procedure we 
assume that the original positions of d non-colliding particles are all differ- 
ent and we consider the limit as these initial positions converge to zero; in 
this way their final distribution is given by a continuous analogue of the for- 
mulas (flOl) and dH]) which give the square of the number of Young tableaux 
of a given shape, with the transition probabilities replaced by the Gaussian 
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kernels. One can easily check that such a derivation of the distribution of 
eigenvalues of a GUEq random matrix follows ([8]) very closely. 

To summarize: our proof of the main result will be based on the observa- 
tion that both Young tableaux and the eigenvalues of matrix- valued Brow- 
nian motions can be interpreted as non-colliding particles and applying the 
determinantal formula of Karlin and McGregor [,KM59,I . 

1 .7. Final remarks. We can see that both the case when d is finite and the 
case considered in Section [L3] corresponding to d = oo are asymptotically 
described by GUE random matrices. It would be very interesting to find a 
direct link between these two cases. 

2. Proof of the main result 
For a function f : M ^ M we define its difference A^J : M ^ M by 

i(y + Ji)-ny) 



By iterating we define A"f for any integer cx > 0. We also define its shift 
: M ^ M by 




n 

Notice that S{^f is well-defined for any integer a. 

Lemma 5. For each n we define a function f : M ^ M which is constant 



on each interval of the form 




(2) fn ^ = 




for each integer k and such that 



ifk is a non-negative integer, 



ifk is a negative integer. 
Then for each integer <x> and y eM. 

(3) lim A^f^(xj) = J=e"*^. 



Furthermore, for each oc > there exists a polynomial Pa such that 
(4) |AX(-y)| <Pjy)e-l^l 

holds true for all n and y. 
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Proof. Before presenting the proof we notice that is a density of a prob- 
ability measure arising as follows: we normalize the Poisson distribution 
with the parameter = in order to have mean and variance 1 and we 



0, 



; therefore 



convolve it with a uniform distribution on the interval 

([3]) states for a = that the suitably rescaled probabilities of the Poisson 
distribution converge to the density of the normal distribution. The case 
oi > 1 shows that this convergence holds true also for differences (respec- 
tively, derivatives). 

The proof of © in the case a = is a straightforward application of the 
Stirling approximation log z! = (z + log z — z + ^^^^ -\- O (z^^ ), namely 

k— — 

for y = such that k is an integer we denote c = Then 



logfn(y) = {c+v^/c+'\^ logo -c -log (c+yVc) ! = 
- (c + yv/^+ 1) log (l + ^) +yv^- + O (C-) 



+ O c-2 



_ _ log 
2 2 

where the above equalities hold true asymptotically for y bounded and c 
oo. 

In order to treat the case a > 1 we observe that the iterated derivative 
on the right-hand side of Q can be calculated by using the following three 
rules: 

-^e^=-ye^; ■j-V = '^; -^i<\>^] = ( A^+^^^-^- 
dy dy dy V^y / dy 

Similarly, the iterated difference on the left-hand side of ([3]) can be calcu- 
lated using the following three rules: 

Anfn = -gnSfn; A^^g^ = 1; A^(ab) = (A^ajb + (Sna)A^b 

where Qn : M ^ M is a function which is constant on each interval of the 



form 



±1 ^+1-1 



for each integer k and such that 



k+1-B 




For each integer (3 we have limn-^oo (Sj^fn) (y ) = 2 and limn-^00 (S^Qn) (y ) 

y therefore each term contributing to the left-hand side of ^ converges to 
its counterpart on the right-hand side of Q, which finishes the proof of dH). 



PERMUTATIONS WITHOUT LONG DECREASING SUBSEQUENCES 



7 




There is a constant Ci < with a property that if y < Ci then the right- 
hand side of ^ is greater than 1 for any value of n. It follows that if 

Vi = -7=- for t G {1 , 2} and -y 1 < < Ci then 

V d 

(6) fn(-yi) <fn('y2)e^'-^^ 

Similarly we find a constant C2 > with a property that if C2 < y 1 < Vz 
then 



(7) U-yz) <fn('yi)e^'-^^ 

For a = inequality dH) holds true for y in a small neighborhood of the in- 
terval [Ci , C2] for Pec being a sufficiently big constant which follows from 
([3]) and compactness argument. Inequality @) holds true outside of the in- 
terval [Ci , C2] by inequalities Q and ([7]). 

The case a > 1 can be proved in an analogous way to the above proof of 
([3]): we show that A^f^^ is a sum of the terms of the form (Sj^^ gn) • ■ ■ (Sj^^ gn) (S|^fr 
and the absolute value of each such a term can be easily bounded by P(y ) e^'^', 
where P is a suitably chosen polynomial. 

□ 



Proof of Theorem\3\ The following discussion is based on the work of Zeil- 
berger nZei83l1 . Every Young tableau T with at most d rows and n boxes can 
be interpreted as a trajectory of d non-colliding particles Xi (t), . . . , Xd(t) 
on the real line as follows. We set 

Xi(t) — d+1 — i+(number of boxes of T in row i which are not bigger than t). 

In other words: the initial positions of the particles are given by (xi(0), . . . ,Xd(0)) = 
(d, d — 1 , . . . , 1 ). In each step one of the particles jumps to the right; the 
number of the particle which jumps in step t is equal to the number of the 
row of the Young diagram T which carries the box with a label t. The condi- 
tion that T is a standard Young tableau is equivalent to xi (t) > ■ • • > x^[t] 
for every value of < t < n. 

Thus the results of Karlin and McGregor nKM59l can be applied and the 
number of standard Young tableaux of the shape Ai , . . . , A^, where |A| — 
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Ai + ■ — h Ad = n, is equal to the determinant 

1 



(8) Na, x^=n\ 



j_ 

Ai! 
1 



(A2-1)! 



(Ai+1)! 
1 



nie 



(5)' 



(Ad-d+1)! {Ad-d+2)! 



:-d+i 



:-d+2 



where 
(9) 



d(d+1 ) 



S/+ifn(yd) AnS/+ifn(-yd) 



(A,+d-l)! 
1 



(A2+d-2)! 



1 

Ad! 



We are interested in a probability distribution on Young diagrams with n 
boxes with the probability of (Ai , . . . , A^) equal to 



(10) 



1 



Cn,d 



(NA,,...,Ad)', 



where Cn,d is the suitably chosen normalizing constant. Clearly, 



Cn.d 



/ X 2T1+-S— ±A£ 



2p2n 



A] ,...,Ad-i 



d-1 



S^^fn(y2) A,S^ifn(y2) 
S/+^fn(yd) A,S/+ifn(yd) 



where the sum runs over Ai , . . . , A^-i such that for A^ = n — (Ai + ■ • ■ + 
Ad-i ) we have that Ai , . . . , A^ is a Young diagram with n boxes. The right- 
hand side can be viewed as a Riemann sum; Lemma[5] shows that the domi- 
nated convergence theorem can be applied (with the dominating function of 
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the form P(-yi, . . . ,'yd)e^^^'^^ l+- +|ydl)^ where P is some polynomial) and 



lim Cr^^a 



(5) 



2n+ 



2c>2n 



e 2 
e 2 



_4 
e 2 



d 

d -Hi 
-jS_e 2 



d 



d-i 



Hi 

2 



Hi 

2 



where the integral runs over (y i , . . . , tf ^-i ] such that for t) ^ = — i H 1" 

ijd-i) we have ui > ■ ■ ■ > ya- 

Since the limit density defines a probability measure, in the limit n — > oo 
the random variables (y i , . . . , y a i ) (please notice that due to the constraint 
y 1 + ■ ■ ■ + U d = the value of y ^ is uniquely determined by y i , . . . , y ^i-i ) 
converge in distribution to the probability measure on the set yi > yi > 
• ■ • >Vd~^ > -(y 1 H h y d-i ) with a density 



C 



Hi 

" 2 



dyi 

d -d 



Hi 

2 



Hi 

2 



d -Hd 

d-Vd 



dv? 

dd- 
dyf' 



dd- 
dyr 



re 



Hi 

2 



Po(yi)e 
Po(y2)e~'^^ 



Hi 

2 

4 



Hi 

2 



Po(Vd)e "2" 



Pi(yi)e 

Pil-yije-^ 



Pil-ydje "2" 



Pd-i(yi)e 2 

_v| 

Pd-i(y2)e 2 



Pd-i(Vd)e "2'' 



for a suitably chosen normalizing constant C^, where -^e^^ = Pklz-je^^ 
for some polynomial pic (related to Hermite polynomials). Since Pk(z.) = 
(— z)'' + (summands of lower degree) the above expression takes a simpler 
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form: 



C 



e 2 
e 2 



_4 



-yi)e 2 
-yije 2 



-ydje" 



^A 



1<i<j<d 



When we set Xt = Vlyi = y ^ {^n,i ~ it becomes clear that the 
limit distribution of (xi , . . . , x^) coincides with the distribution ([T]) of the 
eigenvalues of a GUEq random matrix, which finishes the proof. □ 
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